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ABSTRACT 

0^ . In the innermost region of the advective accretion disk orbiting a black hole of high spin, 

the inertia of heat stored in the accreting gas is comparable to that of the gas rest mass itself. 
, Accounting for this effect, we derive additional terms in the disk structure equations, and show 

. that the heat inertia plays a significant role in the global energy conservation and dynamics of 

accretion in the relativistic advective disks. 
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Introduction 
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If^ • In advective accretion disks (see Narayan 1996 for a review), the released gravitational energy is not 

, radiated away as in the standard model of thin disk (Shakura 1972; Shakura & Sunayev 1973). Instead, it is 

stored in the form of internal heat, being eventually advected into the black hole. Optically thin advective 
disks have been applied to X-ray transients, low luminosity active galactic nuclei, and the Galactic Center 
(Narayan 1996). Accretion can also proceed in the advection dominated regime in the optically thick super- 
Eddington accretion disks (Abramowicz et al. 1988), and the observational significance of these disk models 
2 ' for quasars has been recently discussed (e.g. Szuszkiewicz, Malkan, & Abramowicz 1996). 

-4— » ■ 

^ I Most of the models of the advective disk have been constructed under assumptions similar to those 

used in the standard model: that the flow is stationary and axially symmetric, and that one dimensional 
approximation, with the "integrated" vertical structure, is physically adequate. At the same time, the 
' equations of such a slim disk, differ considerably from the standard ones. In particular, the slim radial 

■ Euler equation accounts for a deviation of gas motion from the Keplerian one, which becomes significant in 

the case of advection dominated disk. A Newtonian version of the slim disk equations has been derived by 
Paczynski and Bisnovatyi-Kogan (1981) long time ago and used by many authors. The relativistic version of 
these equations was derived more recently by Lasota (1994; see also Abramowicz et al. 1996; Abramowicz, 
Lanza, & Percival 1997; Peitz & Appl 1997). Lasota's equations have been integrated numerically to obtain 
relativistic models for the hot disk of low luminosity (Abramowicz et al. 1996; Jaroszyhski & Kurpiewski 
1997; Peitz & Appl 1997; Igumenshchev, Abramowicz, & Novikov 1997). 
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Accreting black holes are expected to acquire a large spin and then the relativistic effects arc especially 
important. A special feature of an accretion disk around the rapidly rotating black hole is that the viscously 
dissipated power becomes comparable to Mc^, M being the accretion rate. As a result, the inertial mass 
associated with heat stored in the advective disk is approximately equal to the rest mass of the accreting gas. 
Lasota's equations, being fully relativistic in all other respects, neglect the inertia of internal heat, and Peitz 
& Appl (1997) and Jaroszyhski & Kurpiewski (1997) have partly included it in their calculations. In this 
note we show that accounting for the heat inertia is necessary to make the disk model consistent with the 
global energy conservation law, and derive the disk structure equations including this effect. In particular, 
an additional term appears in the radial Euler equation. We dub it "heat deceleration" as it describes the 
back reaction of energy release on the radial velocity of accreting matter. 



2. Notation, assumptions and methods 

Our approach and notation are very close to those in Abramowicz et al. (1996), except that we do 
not neglect the contribution of internal energy and pressure to the inertial mass of the flow. We denote 
by rg = 2GMI<? the gravitational radius of the black hole with the mass M, and by a = J/Mc the Kerr 
parameter connected to the black hole angular momentum J. We use the Boyer-Lindquist coordinates 

a;* = {t,r,6,(p), and the signature ( h ++). The metric tensor gij is given e.g. by Misner, Thome, & 

Wheeler (1973). The four-velocity of the accreting gas has components = (u',u''\u^,vy) in the Boyer- 
Lindquist coordinates. The disk is assumed to lie at the equatorial plane of the Kerr geometry. In the slim 
approximation, is neglected and all the metric and connection coefficients are evaluated at the equatorial 
plane. 

The gas Lorentz factor measured in the frame of local observers with zero angular momentum is con- 
nected to u* by 7 = u*(— g**)^^/^. The angular velocity of the gas rotation is defined as O = u'^/u*. It is 
not equal in general to the Keplerian angular velocity, 

fi±=± E 

^ r(2r/rg)V2±a" 

(Plus and minus signs in this formula correspond to co-rotating and counter-rotating orbits). The r- 
component of four-acceleration, Ur, is given by (Abramowicz et al., 1996) 

ar = ~ {uru-) + i^5"7^ - ^i) {n - n-) . 

The dynamical equations are derived from the conservation laws ViT'^f. = 0, where V, is the covariant 
derivative in the Kerr metric, and T'^. is the stress-energy tensor of the viscous gas flow (Misner et al. 1973), 

T% = {p + p)u'uk +p5\ - 2cT]a\ + ^{u'qk + Ukq'), 

where p and p are the total energy density and pressure in the comoving frame, (T*^, is the shear tensor, and 
is the energy flux vector assumed to be directed vertically in the disk. The dynamic viscosity 77 is related 
to the kinematic viscosity u hy rj = iy{p + p)/c^. (Note that this differs from the previously used relation 
that includes only the rest mass of the gas). 

The flnal equations are written for the vertically integrated thermodynamic quantities: surface rest mass 
density S, surface energy density U (that includes both the rest mass energy Sc^, and internal energy 11), 
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and the vertically integrated pressure P. F+ denotes the surface rate of viscous heating, and F~ denotes 
the radiation flux radiated from both faces of the disk. All the thermodynamic quantities, and both fluxes 
F~ , are measured in the comoving frame. The dimensionless specific enthalpy is defined as, 

U+P U+P 

Neglecting the internal heat contribution to the inertia of the flow is the same as assuming that = 1. In 
this paper wo do not assume that, but keep in all equations fi > 1. Consistently, we keep the viscous term 
in the radial Euler equation (see Section 4). 

We skip all computational details, because our derivation goes along the same standard lines explained 
e.g. by Page & Thorne (1974), and used recently by Lasota (1994) and all the other above quoted subsequent 
authors. In the next Sections we give only the final results. 



3. Conservation of energy and angular momentum 

The conservation of energy and angular momentum are expressed by equations 



d_ 
dr 

d_ 
dr 



rut, 



Mu. 



ru 



Here M is the accretion rate related to and S by the barion conservation law, 



-M, 



(1) 



(2) 



(3) 



and a^f, a^^ are the shear tensor components, 



dQ 
dr ' 



These equations differ from those with neglected inertia of heat (see Abramowicz et al. 1996) just by the 
presence of the factor /i > 1. They also differs from the equations in Peitz & Appl (1997) and Jaroszynski & 
Kurpiewski (1997) who take into account the contribution of heat to the specific orbital energy and angular 
momentum of the flow but neglect the similar term in the relation between the dynamic and kinematic 
viscosity. 

In the standard thin disk /it w 1 with high accuracy, and the rotation is very close to Kcplerian. Then 
the equations (1,2) become the same as those derived by Page & Thorne (1974). For the advective disk, 
however, the factor fi should not be neglected. This can be easily seen from the global energy conservation 
law. Integrating equation (1) from the transonic inner edge of the disk, rj„, to infinity, and neglecting the 
viscous stress at ri„, we get 



L = 



2-K 
c 



utF~rdr = Mc^ 1 + M 



(4) 
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It follows that the radiative efBciency of the disk equals 

e = 1 - M»n ^1 - , (5) 

where bin = + ul"c is the specific binding energy at the inner edge. In the advective disks, most of the 
dissipated binding energy is stored as internal heat and the radiative losses are negligible, e — > 0. ^^From 
this we conclude that ^ (1 — bin/c^)^^- In the particular case of the extremely rotating black hole, 
assuming e.g. that bin/(? is close to the corresponding standard value (1 — l/\/3) ~ 0.42, one has jjLm = VS. 
Note that the usual assumption ^ = 1 follows that bin = for any advective disk. In fact, the position and 
binding energy of the inner edge can be found only by the integration of the disk structure equations, and 
it is the difference (/x — 1) that adjusts to keep e « 0. In this respect the heat inertia is essential in the case 
of non-rotating black hole as well as extremely rotating. 

It would be instructive to compare the angular momentum equation (2) with its commonly used standard 
version (Novikov & Thorne 1973; Page & Thorne 1974). While ^ = 1 can be safely assumed in the standard 
disk, the radiative losses of angular momentum represented by the right hand side of equation (2) become 
important when the accreting black hole is rapidly rotating {e.g. Lamb 1996). Contrary to this situation, 
in the advection dominated disk the radiative losses are always small, but the deviation of fi from unity 
becomes significant. In a sense, the angular momentum which would be radiated away by the standard disk 
remains now in place, being carried by the increased inertia of internal energy. 



4. Radial motion and viscous heating 

To clarify the role of heat inertia in the equation of radial motion we derive this equation in two ways. 
Firstly, we perform vertical integration of the r-component of the equation VjT'j, = 0, to get 



U Ur 



dP F~ 

ar{U + P) + — +Ur =0, (6) 

dr c 



where ^ « 1 is a numerical factor accounting for non-homogeneity of the disk in vertical direction. In this 
equation we have taken into account that the divergency of the vertical energy flux becomes equal to F~ 

after the vertical integration. 

Secondly, we project the equation ViT^^ — onto the hypersurface orthogonal to to get the relativistic 
Eurler equation (c./. Lightman et al. 1979, problem 5.31). Then the flux term vanishes, but the viscous 
term arises in the radial equation, and we get after the vertical integration 

dP F+ 

ar{U + P) + —{l+vrUr) + Ur=0, (7) 

dr c 



where 

F+ = 2^Em a'c\ = \g^^ g^^ (-5") 7' V • (8) 



Keeping terms in the equations (6,7) allows to get the first law of thermodynamics for a vertically 
integrated accretion disk as a consequence of these two equations. 
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Finally, substituting a^- to (7), we write down the equation of radial motion with u'^Ur <C 1, 

5 i = -1 (" - "a (" - "i) - - Sir- m 

This equation differs from the radial equation in Abramowicz et al. (1996) by the presence of the factor 
fi, and by the additional term proportional to . Without this term, equation (10) would have a simple 
meaning: radial acceleration is a combined result of the deviation from Keplerian rotation and the radial 
pressure gradient. The viscous term oc represents the heat deceleration effect which can be interpreted in 
the following way. The mass of the accreting gas measured in its local rest frame increases due to the stored 
heat. Thus, from the local observer point of view, the dissipation of orbital binding energy is an external 
source of mass-energy, and the matter inflow proceeds like motion of a body with changing mass. In this 
case, in addition to acting forces there is a contribution to acceleration due to the change of mass. The mass 
of the accreting gas increases with rate proportional to and this tends to decelerate the matter inflow. 
The "deceleration by heating" is a purely relativistic effect as it accounts for the mass-energy relation. It 
is proportional to the ratio of the released power to Mc^ and can be essential near a rapidly rotating black 
hole. Then the additional term must influence the structure of the advective disk in the innermost region. 
There it becomes of the same order as the left hand side of the equation (10). 

Note that the deceleration term is independent of the radiative losses, F~ , and remains the same even 
when F~ ~ In this case the bulk of the dissipated energy is radiated away. However, this energy first 
appears as heat decelerating the radial inflow. Only then it is radiated away. The net momentum taken 
away by the radiation flux vanishes in the rest frame of the accreting gas, and does not contribute to the 
velocity change. 
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